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Abstract 



We present a model to find analytically the electronic states in self-assembled 
quantum dots with a truncated spherical cap ('lens') geometry. A conformal 
analytical image is designed to map the quantum dot boundary into a dot 
with semi-spherical shape. The Hamiltonian for a carrier confined in the 
quantum lens is correspondingly mapped into an equivalent operator and its 
eigenvalues and eigenfunctions for the corresponding Dirichlet problem are 
analyzed. A modified Rayleigh-Schrodinger perturbation theory is presented 
to obtain analytical expressions for the energy levels and wavefunctions as 
a function of the spherical cap height b and radius a of the circular cross 
section. Calculations for a hard wall confinement potential are presented, and 
the effect of decreasing symmetry on the energy values and eigenfunctions of 
the lens-shape quantum dot is studied. As the degeneracies of a semi-circular 
geometry are broken for b ^ a, our perturbation approach allows tracking 
of the split states. Energy states and electronic wavefunctions with m = 
present the most pronounced influence on the reduction of the lens height. 
The analytical expressions presented here can be used to better parameterize 
the states in realistic self-assembled quantum dots. 

PACS99: 73.61.-r, 73.20.Dx, 03.65.Ge, 78.30.Fs 
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I. INTRODUCTION 



Quantum dots obtained by interrupted growth in strained semiconductor interfaces are 
currently under intense study by many experimental and theoretical groups.§ These 'self- 
assembled' quantum dots are mostly dislocation free, coherent islands of deposited material 
on the surface of a different semiconductor. The lattice mismatch from one semiconductor 
to the other forces the segregation of material whenever the epitaxial growth exceeds a crit- 
ical layer thickness, resulting in the growth of these so-called Stranski-Krastanov islandsB 
Deposition of material past a critical thickness, which depends on the two materials used, 
results in large arrays of small islands with a rather narrow size distribution (with size vari- 
ations well under 10%), and typically arranged randomly on the plane (although avoiding 
overlapping islands, for the most part).0 More recently, some groups are working at pro- 
ducing in-plane ordering of islands following different approaches, including 'nucleation site 
engmeermg ■ to favor certain locations for =df-a=sembled dot growth! These self-assembled 
quantum dots have, for the most part, a large area-to- height aspect ratio. In the case of 
InAs islands grown on a GaAs surface, the in-plane diameters are typically less than 30 nm, 
while their heights are below 10 nm.S Although not without some controversy, it is generally 
believed that the InAs quantum dots on GaAs are shaped like a 'lens', characterized by a 
spherical cap shape with circular cross section. Upon optical or other mechanisms of carrier 
injection, both electrons or holes have bound well-confined states inside these dots. 

In fact, optical experiments on self-assembled dots demonstrate that these structures 
provide strong carrier confinement, as decreasing dot size produces strong blue-shift of the 
extremely-narrow luminescence peaks in isolated dots.i^ Confinement effects have also been 
shown to appear in magneto-capacitance and infrared absorption experiments by several 
groups.llytJ Clearevidence of electronic shell states and their different degeneracies has been 
reported recently,lij and novel applications such as storage of photoluminescence signals was 
demonstrated recently in gate-activated 'optical memories. 'tJ 

Simplifying the symmetry of the lens as a two-dimensional harmonic oscillator has been 
suggested to characterize the level structure for charge carriers in the dots, and this has 
proven useful in the interpretation of experiments.0 Nevertheless, as we will show here, a 
lens geometry has quite a different level structure and wavefunctions with subtle symmetries 
which might be seen in experiments. For example, we find that as the height of the cap or 
lens decreases, there is a larger shift of the wavefunction towards the plane of the lens (in 
comparison with the situation in a semi-spherical geometry) . This shift becomes stronger for 
'flatter' lenses, and may even give rise to deconfinement of the state towards the substrate, 
for finite confinement potential, changing significantly the oscillator strength of electron- 
hole transitions, for example.0 However, even before that occurs, the smaller height lens 
geometries exhibit a different level structure than the harmonic oscillator, as we show below. 
This structure may give rise to different Pauli blocking effects and transition rules, topics 
which are the subject of interest in recent experiments.ll3 

In this paper we rigorously show that a complete set of wavefunctions can be generated 
with the correct lens symmetry to describe the physical properties given by the differential 
equation of interest. For the sake of simplicity, we have applied our calculation here to the 
single-particle Schrodinger equation. The incorporation of non-parabolic band dispersion or 
many-particle interactions is straightforward (if only a bit cumbersome). Notice that while 
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our calculation provides an analytical approach to the understanding of electronic states 
in these structures, it also provides an interesting example of a generalization of perturba- 
tion theory for the conformally-mapped differential operator arising from the Schrodinger 
equation in the material. Since the conformal transformation used to map the lens into the 
semi-spherical geometry is non-trivial, the resulting differential equation in the transformed 
space reflects that complexity. Fortunately, our perturbation approach is robust for the 
identified small parameter of the problem, as we will show. 

The remainder of the paper is structured as follows: in section II we introduce the 
problem of the spherical cap geometry, its semi-spherical limit, and the conformal map 
which connects the two shapes. In that section, we also describe the perturbation approach 
needed to carry out the solution of the appropriately mapped Schrodinger operator, and 
explore the orthonormalization conditions of the basis used in the description of the general 
problem. In section III we present some examples of the eigenstates for varying lens size, and 
analyze their angular and radial distribution functions. Section IV presents some discussions 
and conclusions, while the appendix contains details of the expansion described in section 
II. 



II. DIRICHLET PROBLEM FOR A QUANTUM LENS: CONFORMAL MAPPING 

The shape of the quantum lens in real r-space is shown in Fig. 1(a). R^{a,h) denotes 
the domain of the lens with boundary in r-space given by a spherical cap of height h and 
circular cross section with radius a. Assuming a carrier Hamiltonian model with isotropic 
band and effective mass m*, the eigenvalue problem for a particle confined in the lens is 
described by the operator 

^ = -^V^ ' r-Gi?3(a,&), (1) 

obeying the Dirichlet boundary condition ip = Q for r G L^{a, b). L^{a, h) is the boundary of 
the R^{a, b) domain. The operator (|I]) presents axial symmetry and all functions defined on 
i?3(a, b) have the property = (p{ip + 27i), where (p is the axial angle. Hence, the solution 
of (|ip can be written as 

^{r) = f{p)^= ; m = 0,±l,±2,..., (2) 
V 27r 



where p is a 2D vector. Correspondingly, the operator (|l|) is transformed into an eigenvalue 
problem for /(p), 

'''"'/(p)=0 ; pGi?2(a,6), (3) 

where = 2m*E/h^, E is the energy of the eigenfunction /(p), 9 is the polar angle, 
and R2{a,b) is the 2D domain with boundary L2{a,b) in the p-space shown in Fig. 1(b). 
The Hilbert space where the operator (0) is defined corresponds to the set of functions 
/(p) G i?2(a,fe) with boundary condition /(p) = on L2{a,b). 
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A. Semi-spherical quantum dot 



For the case of b = a, the lens shape of Fig. 1(a) has semi-spherical symmetry. Hence, 
the Dirichlet problem in (|^) reduces to the conditions f{a,9) = 0, and /(p, 7r/2) = 0. The 
set {fi} of eigenf unctions of on the R2 domain forms an ortho normal basis with functions 
which in polar coordinates are given by products of the associate Legendre polynomials and 
Bessel functions, 

= (4) 

with / = 0, 1,2, and the condition — I < m < I. The normalization constants Ni^m and 
Nb are as usual given by 



N, 



m 



\2(+l ((-|m|)! • ° V2 



where is the derivative of the Bessel function J^, and //^^^ is the n-th zero, Jp+i(/i^''^) = 

0. The eigenvalues are given by En^i = fi^ (/^n'*) /(2m*a^), and the boundary condition 
f^^\p,9 = 7r/2) = restricts the values of the quantum numbers / and m to fulfill the 
condition \l — m\ = odd. According to this condition, the degeneracy of states f^] for a 
given energy En,i is equal to / and the ground state corresponds to / = 1, m = 0, and n = 1. 



B. Quantum lens 



The quantum dot with lens shape corresponds to the more general case when b < a. 
Here, we need to fulfill Eq. with the Dirichlet condition over the boundary L2{a, b). The 
wavefunctions f^] given by (^) are not solution for the general case, because the problem has 
no longer the semi-circular symmetry in 6. The energy number n and angular momentum I 
are clearly no longer good quantum numbers when b ^ a and the m-degeneracy is broken. 
To obtain an analytical solution of the problem it is convenient to make a conformal 
mapping to the circular cap with domain R2{a,a) and boundary L2{a,a). The mapping 
enables us to solve the Dirichlet problem for the operator given by in a Hilbert space 
where an orthonormal basis {fi} is known. We transform the quantum lens domain and 
boundary into a dot with semi-spherical shape, so that the circular cap defined by the domain 
Z = X — iz G R2{a,b), transforms into the semi-circular domain W = u — iv G R2{a,a). 
This is accomplished by the transformation (see Fig. 1(b)) 

~ 1 I ( a-Z \^ ~ ^ ' " ~ arctan(6/a) ' 
\a + Z) 

In the complex plane W, we have the parameter equations: u = psinO, v = pcos9, < p < 
a, < 6 < n /2. The eigenvalue problem (^) is thus transformed by this conformal mapping 
into the problem 
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V?„,,)F(u, v) + Ja{u, v)[e- ' F{u, v) = 0; {u, v) G i?2(a, a), (7) 

with the boundary condition, 

^(«'^)l(M,i;)eL2(a,a)=0- (») 
dZ ^ 



The functions J7a(u, f ) = (the Jacobian) and Jaj are given in the Appendix, and 

are mathematical objects which contain the information of the lens geometry, where the 
subscript a is given in (^). It should be noted that a > 1, since b < a, and for a = 1 the 
Jacobian JT^ reduces to 1, while Jaj^a reduces to Xju^ = 1/ (psin6')^. 

The Hilbert space on which the operator (|^) is defined must fulfill the Dirichlet boundary 
conditions indicated in (P). A set of functions that fulfill these boundary conditions are the 
functions f^j given in (Q), and represent a complete set of orthonormal eigenfunctions for 
the operator in The solution F{u,v) for a given m can be expanded in term of the set 
{fnd} such that 



F = T.Cn,ifSiP)^ (9) 



n,l 

where p= {p,6) is here the parameterization of (u, f), and the functions f^] are restricted 
to the condition |Z — m| = odd. The coefficients Cn,i have to be determined to satisfy the full 
operator; a perturbation procedure to accomplish this is described in section |11D| below. 



C. Orthogonality and completeness 

Equation can be cast in operator form as 

KF = J^k^F, 



(10) 



where K involves the Laplace operator, — V(„„), and the term {m? — 1/ A)Ja/ X^{u,v). 
Equation (p!0D is an eigenvalue problem for the dimensionless energies klflm), where is a 
generic label for the different eigenstates Fn of (pi][). Let us now suppose that A^ and M 
correspond to different eigenvalues of ([T0|). From the above equation it follows that 



f;jKf^ 



FnK F 



M 



^^m)FmFn . 



(11) 



Making the integration over p{u,v) in the R2{a,a) domain, we obtain that 



Ir 



IR2(a,a) 

Integration by parts gives us 



Fm^Iu,v)Fn + f^v'mf;, d'p = (4 - k 



{u,v) 



''MJ 



R2{a,a) 



Fl^F^Jad^p. (12) 



(4 - kll) I JcFIiFm d^ = [FnWFI, - FIjWFn) |(„,.)gL,(a,a) 

J R2{a,a) 



(13) 
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which due to the boundary condition (H), it is reduced to 



R2{a,a) 



M 



FNd^P 



0. 



(14) 



For 7^ M, condition (T^) represents the orthogonahty property of the eigenfunction set 
{Fn}, where jTo is clearly the weighting factor of the eigenproblem (0). Moreover, the 
operator K is Hermitian, ensuring that the solution of the present problem is described by 
means of a complete orthonormal basis of eigenf unctions {-F/v} obeying the condition 



R2{a,a) 



JaFljFN d^p 
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N,M 



(15) 



D. Perturbation theory 

The coefficients C„ ; in (|) and the eigenvalues fc^ can be obtained by perturbation theory 
if 6 ~ a (a — > 1). In this case, the lens cap represents a perturbation from the semi-spherical 
geometry. In other words, the operator (0) can be rewritten in the form 

{Ho + Hj,)F{u,v) = Q, (16) 

with 



H,{u,v) = e{Mu,v) - l)-(m^ - 1/4) - - . (18) 

The operator Hp vanishes when a ^ 1 {b ^ a) and it can be considered as a small 
perturbation operator, when the height b is close to the radius a. The set {fn]} given by Eq. 
(I) are the eigenfunctions of the Hamiltonian Ho in the W-space and form an orthonormal 
basis on the R2{a,a) domain. In order to find the solution of (|1^) as a function of the 
ratio b/a, we will develop a modified Rayleigh-Schrodinger perturbation theory. We note 
that the perturbation Hamiltonian Hp depends on the eigenvalue k'^, and as such requires a 
somewhat different approach. Substituting Eq. in ([T6|) we obtain 

\k'-kl) + {n,l\Hp{e)\n,l)]Cr,,i + 

+ {n,l\Hpie)\n',l')C^,,i, = 0. (19) 

n',l'jtn,l 

The states of Ho are degenerate on the quantum number m. Nevertheless, according to (^, 

(n, l,m \ Hp \ n', I', m') = {n,l \ Hp \ n', I') 6m, m', (20) 

so that the matrix elements are diagonal on m and we can develop a perturbation theory 
in the absence of degeneracy. We represent the coefficients C„ ; and its eigenvalues k^ in a 
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power series of the small parameter \ = a ^ — 1 (which arises naturally from the expressions 
in the Appendix). We obtain up to first order in A an expression for the wavefunctions given 
by 



2 \ ' 



dX 



A=0 



n,l) fi']{p,9) 



n', I' n,l 



kl{n' ,1') - klM 



X 



n'A' 



dX 



[m^ - 1/4) — 



d f Ja 



A=0 



dX \ 



A=0 



n, / 



We find for the eigenvalues up to second order, that 

k'^{N, m) = kl{n, I) + X kl{N, m) + kl{N, m), 



where 

kj{N,m) 

kliN,m) 



n, I 



kl{n,l) 



E 

n', V ^ n^l 



\ dX 
n', I 



[m' - 1/4) — 



d ( J. 



A=0 

KM 



dX y^a^ J x=0 



n, / 



(21) 



(22) 



(23) 



~5X 



A=0 



a / A=0 



n, I 



-- ( nj 



kl{N,m) 



dX 



-{m' - 1/4) 



A=0 V 
dX^ V'^a /A=0 



kl{n\l') - klM 



dX^ 



n, I 



A=0 



(24) 



The different factors included in the geometric perturbation, 



dX 



A=0 



dX^ 



A=0 



dX V XJ 



" / A=0 



5A2 \ X, 



(25) 



" / A=0 



are also given in the Appendix. In the framework of the infinite confinement model, the 
parameter dependence of Eqs. (pTD-(p^) is known, since the matrix elements play the role of 
constants and need to be evaluated only once. Note also that the above expressions are not 
the same as those found in a typical perturbation theory, as the difference arising in the k"^ 
term depends on the perturbation Hamiltonian Hp, which itself depends on the parameter 
A. 



III. THE EIGENVALUES AND WAVEFUNCTIONS 

Figure 2(a) shows the first 13 energy levels in units of Eo = h"^ / {2m* a?) for a quantum 
lens calculated by perturbation theory up to the second order in A = a"^ — 1, as a function of 
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the ratio h/a. The different eigenvalue curves are labeled by the quantum numbers {N,m). 
The semi-sphere case {b/a = 1) is indicated by the limiting value on the right vertical axis 
in each panel. One can see the breaking of degeneracy in the quantum number m, and the 
strong deviation from the semi-spherical the ratio b/a decreases. The lower levels 

exhibit a weaker dependence on the decreasing b/a ratio, while the upper levels are strongly 
deviated from the semi-spherical case. In Fig. 2(b) the first 5 energy levels calculated up 
to first (dotted lines) and second order (solid lines) perturbation on the parameter A are 
compared in the range 0.4 < b/a < 1. It can be seen that a strong deviation is present for 
the higher excited levels {N = 3 and 4), while for N = 1 and 2 the obtained results using 
Eq. ( [23| ) (first order perturbation theory) give a deviation smaller than 1% in comparison 
with those using Eq. ( p4|) (second order perturbation results). 

The radial and angular probability density function in a given state A^, m are defined as: 



t/2 a 



(26) 



The functions P{p) and P{0) are obtained up to first order perturbation theory on A ac- 
cording to the equations 



PN,m{p) = P 



' Ji+i/2ifJ'i^Jp/a 



N, 



B 



and 




A=0 



nA) + 



P e (0,a) 



(27) 



P^,^(^) = sin^'^'""'^'°" 



+ 



\ aA 




A=0 



l + \y-{n,l 
'Ji+i/2{pll^p/a 



B 



dX 



A=0 



nJ) + 



; (0,7r/2). 



(28) 



Figure 3(a) and (b) show the radial and angular probabilities given by Eqs. ( P7[ ) and 
respectively. In both cases the levels considered are {N,m) = (1,0), (4,0), (8,0), and 
(8,2). Solid lines represent the semi-spherical case {b/a = 1), while the quantum cap with 
b/a = 0.509 is shown by dotted lines, and illustrates the departure from the semi-spherical 
case for decreasing b/a ratio. In the case of the radial probability, the deviation observed 
when the ratio b/a decreases is relatively small, in comparison with the semi-spherical cap. 
In contrast, the angular probability shows a rather strong deviation as a function of b/a, 
as the maximum probability is shifted towards 6 = tt/2; that is, the carrier is located more 
towards the plane f = 0, as the cap height decreases (see Fig. 1(b)). For all states, the 
maximum of P{p) is smoothly shifted towards p = 0, except for level (8,2). The different 
behavior observed for the radial and angular probability densities can be seen as arising from 
the fact that the geometry of the quantum dot is essentially decreasing in radius and not 
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in angle, as h/a decreases. The radial probability is obtained by integration along the angle 
in all directions (0 < 9 < vr/2), taking into account the angular contribution for a given p 
and certain geometry. On the other hand, the angular probability is calculated by taking 
the integration along the radius where the change of geometry is more important. Hence, 
one can say that the angular probability density P{9) contains more information about the 
changing cap geometry than the function P{p), as function of b/a. 

Finally, a test for the viability to obtain the wavefunctions and energies by the per- 
turbation method, is given by the ratio between the matrix element of the perturbation 
Hamiltonian {'n,l\Hp{J'a,m)\n,l) with respect to the non-perturbed dimensionless energy 



The parameter A^v,™. was calculated for the levels of Fig. 2 and is shown in Fig. 4 as a 
function of b/a. A necessary (but not sufficient) condition for the perturbation theory to be 
valid is that AN,m. must be less than the unity, and this criterion is fulfilled for the range 
0.4 < b/a < 1. Another more restrictive condition for the applicability of this method is 
that {n,l\Hp{J'a,m)\n' ,1') < \k'^{n,l) — k'^{n',l')\, saying that the differences between non- 
perturbed dimensionless energy states need to be larger than the matrix elements of the 
perturbation. A more complete, non-perturbation, solution of Eq. (|^ is needed to completely 
assess the validity of the perturbation approach. However, one can be confident that for 
small A values, the eigenstates and eigenfunctions found represent an accurate solution to 
the problem, as the perturbation and the small parameter is well defined, and the procedure 
robust. 



We have presented a formal and systematic conformal analytical map model to describe 
quantum dots with lens geometry and circular cross section arising in the growth of low 
dimensional semiconductor systems. The reported transformation can be extended straight- 
forward to different physical-mathematical models as electronic states and phonon modes 
fulfilling a certain differential equation. The proposed conformal image maps the lens bound- 
ary into a dot with semi-spherical shape, allowing one to obtain a complete set of orthonor- 
mal functions to characterize the physical problem keeping the full lens symmetry. We have 
applied the formalism to the eigenvalue and eigenfunction of the Schrodinger problem in 
a spherical cap geometry. The conformal mapping of the equation allows a modified but 
well defined Rayleigh- Schrodinger perturbation approach, where the cap height to in-plane 
radius is used to define the small parameter of the theory. 

We find that the wavefunctions are strongly shifted towards the flat face, as the height of 
the lens decreases, while the radial dependence is not affected as much. This change in the 
wavefunctions is interesting on its own, as it reflects the changes produced by the appropriate 
operator after the conformal mapping. Moreover, these changes may have important con- 
sequences for the different electronic and optical properties of self-assembled quantum dots. 
We are currently studying those properties and will report our findings in the future. The 






(29) 



IV. CONCLUSIONS 



9 



reported energy dependence on cap height to in-plane radius and semiconductors parame- 
ters can be useful to characterize the geometrical dimensions of these novel semiconductor 
nanostructures. 
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APPENDIX 

The Jacobian of the transformation yV{Z) is given by 

16 (l/«)2 



JJr, 



Jll-l/a jV- ^ ji/" ^ 2/2V2« cos(0/ 



a] 



2 ' 



and the term Jaj'^a can be cast as 



16 {Ijaf 



2 • 



In the above equations we have defined 

r = p/a ; /± = 1 + ± 2r sin ^ ; R = f+f-, 

and 

'arctgf^^") ;r<l 



tt/2 



r = 1 . 



From Eqs. (KTl) and (^), it follows that 



X=i(r,^) = l 



JJr, 



2 sin^ e ' 



as one would expect. 

Finally, the geometric perturbation factors in Eq. p5|) are given by 



'dMr,9y 



d f JJr, 



A=0 



2-rsm9ln{U/f^) + 2(j)rcose. 

4rsin0- (l + r2)ln(/+//_) 
2r3 sin^ 9 ' 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 
(A7) 



10 



A=0 



1 (dJa 

2 I dX 



A=0 



dX 



A=0 



(l + r^) 



In^ 



+02(l-r2) + 21ii2i?, 



A=0 



(4rsing)^ - 16r sin g(l + r^) ln(/+//_) 
8r4 sin^ ^ 

, (2(l + rY + i?)ln'(/+//-) 
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FIGURES 



FIG. 1. (a) Lens cap of height b and radius a. (b) Conformal mapping from Z to W. 

FIG. 2. Energy levels E]\f,jn{b/a), labeled by {N,m), for a quantum lens as a function of 
the ratio b/a. The energies are given in units of Eq = 2m*E/h^. (a) The first 13 energy levels 
calculated up to second order in perturbation theory, (b) Comparison between results calculated 
up to first (dotted lines) and second order (solid) perturbation theory for the first five levels. 

FIG. 3. (a) The radial probability density in a quantum lens PN,m{p), in units of a~^, for 
different electronic states (iV, m), and as function of the dimensionless coordinate r = p/a. (b) 
Angular probability density -P/v,m(^) as function of the angle 6. Two values of the ratio b/a are 
considered: b/a = 1 (solid lines), and b/a = 0.509 (dotted lines). The calculations were based on 
the perturbation theory described by Eqs. (p7|) and (p^)). Different states {N,m) = (1,0), (4,0), 
(8,0), and (8,2) are indicated. 

FIG. 4. Ratio between the matrix element A and the energy k'^{n,l) with respect to the ratio 
b/a of the quantum lens. For all b/a between 0.4 and 1 the parameter A is small. 
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